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Broken time reversal symmetry (BTRS) in d + id' as well as in d + is superconductors is studied 
and is shown to yield current carrying surface states. We evaluate the temperature and thickness 
dependence of the resulting spontaneous magnetization and show a marked difference between weak 
and strong BTRS. We also derive the Hall conductance which vanishes at zero wavevector q and 
finite frequency u, however at finite q, u) it has an unusual structure. The chirality of the surface 
states leads to quantum Hall effects for spin and heat transport in d + id' superconductors. 
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I. INTRODUCTION 

Recent data on the high T c superconductor Y Ba,2Cuj ) O x (YBCO) has supported the presence of broken time 
reversal symmetry (BTRS) 12 3 . A sensitive probe of BTRS are Andreev surface states. For a d wave with time 
reversal symmetry bound states at zero energy are expected for a surface parallel to the nodes (i.e. a (110) surface 
in YBCO). When BTRS is present, by either a complex order parameter or by an external magnetic field, the bound 
states shift to a finite energy. Indeed tunnelling data usually shows a zero bias peak which splits in an applied field; 
the splitting is nonlinear in the magnetic field, indicating a proximity to a BTRS stated. In fact, in some samples 
tunnelling data shows a splitting even without an external field 1,2 , consistent with BTRS; the splitting increases with 
increasing overdoping 2 ^, 

Further support for a spontaneous BTRS state are spontaneous magnetization data as observed in YBCO^, setting 
in abruptly at T c and being almost temperature (T) independent below T c . The phenomenon has been attributed to 
either a d x 2_ y 2 + id xy state {d + id') or to formation of it junctions. No microscopic reason was given, however, for 
the spontaneous magnetization being independent of both T and of film thickness 3 -. 

It has been shown theoretically that BTRS can occur locally in a d x 2_ y 2 superconductor near certain surfacesSiS^ 
leading to either d + id' or d + is states with surface currents. The onset of such BTRS is expected to be below T c 
and therefore does not correspond to the spontaneous magnetization data™. We note that in response to an external 
magnetic field the surface states are paramagnetic and compete with Meissner currents. This effect has been proposed 
to account for a minimum in the magnetic penetration length 10 . In fact, it was proposed that this paramagnetic effect 
leads to spontaneous currents and BTRS in a pure d x i_ y i state 1112 . The onset of this BTRS is much below and 
therefore does not correspond to the data-. 

Of further theoretical interest is the relation of the BTRS state to quantum Hall systems with a variety of 
Hall effectsi 3 ^^^^^. In particular a finite charge hall conductance has been suggested^ 3 *, though this has been 
questioned^. 

In the present work we expand our earlier worki& and study variety of phenomena related to surface currents. In 
section II we show that bulk d + id' state has surface states with finite surface currents; a similar situation was found 
for the bulk p wave stated. We also consider a d + is state which has surface currents only on the (110) surface. In 
section III we evaluate the spontaneous magnetization and show that for d + id' it is dominated by (100) surfaces; 
for thin films it increases with the ratio A'/A (A and A' are the amplitudes of d x i^ y i and d xy , respectively) while 
for thick films it has a maximum at A/£' ~ 1 where £' = Up/A' with vp the Fermi velocity and A is the penetration 
length, i.e. at A'/A ss 0.01 for YBCO. We show that for weak BTRS, A/£' < 1, the spontaneous magnetization is 
T and thickness independent, while for strong BTRS thickness and T dependence may occur. For the sample of Ref 
H we estimate A'/A « 10~ 4 , i.e. weak BTRS. In section IV we consider a surface approach for the quantum Hall 
effect, showing quantization for spin and thermal Hall conductances for the d + id' state. We also derive in section 
V the effective action in the bulk and identify the Hall coefficient which has an unusual wavevector and frequency 
dependence. 

II. SURFACE STATES 

We present here the Bogoliubov de-Gennes (BdG) equations for quasiparticles in a bulk d + id' or d + is states 
in presence of a boundary and study the resulting surface states. We consider first a d + id' state where the order 
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parameter is 

A(p x ,Py) = A(pl-p 2 y )/k 2 F +iA'p x Py/k 2 F (1) 

where p = —ihV is the momentum operator and kp is the Fermi momentum. The quasiparticles are represented by 
an electron-hole Nambu spinor 

/* T (r)\ 

m = + ( 2 ) 

and are described by the following mean field Hamiltonian (see appendix A) 



7-/ = ^-[(iV + ^T 3 A(r)) 2 -4]T 3 + ( | ( -"u- lr . ct) 




where m is the electron mass and Tj are the Pauli matrices. We assume here that |V0| <C kp so that the issue 
of gauge invariance in the interaction term can be avoided (appendix A). Rotating by the unitary transformation 
*(r) — > exp[iT 3 9(r)/2]^>(r) yields 



H = J-(-V 2 -4)T3 + -^-p-(Ve--A(r)) + j 1 (4) 

2m 2m c 



A(p x ,p y )' 

,A*(p x , P y) 



where A is kept to first order. 

We consider a vacuum-superconductor boundary at x=0, and assume for now that A, A' are constants at x > 
and vanish at x < 0. For A > A' this corresponds to a (100) surface; to describe a (110) surface A and A' need 
to be interchanged. The spinor wavefunctions for the up and down component of Eq. l(2|l. respectively, u(r) = 
uexp[ifx + ik y y] and v(r) — vexp[ifx + ik y y] with eigenvalues e satisfy the BdG equations 

if 2 -k F + kl-2mi)u + 2mA(f,k y )v = 
(-/ 2 + k 2 p-k 2 y - 2m~e)v + 2mA*(/, k v )u = (5) 

where e = e + (e/mc)k y Ay(x), A has only an A y component consistent with a current in the y direction and V0 = 0. 
This Doppler shift assumes that A y {x) is slowly varying on the scale kp 1 so that a local eigenevalue e can be defined. 

Define k — +^Jk 2 F — ky, then / has two surface solutions with 3/ > 

. m 



h= k + tj X /\A(k,ky)\ 2 



h= -k + ijy/\A(-k,ky)\*-i? (6) 

where the replacement A(f,k y ) A(±k,k y ) is valid for |A|,e <C k 2 F /2m. The eigenvectors are 



.y/\A(k,k y )\2-?+ii 

Vl = ~ l WJy) ^ 



. y/\A(-k,ky)\*-e>-ie 
V2 = 1 AHUy - 2 ( 7 ) 

We assume specular reflection which preserves k y but mixes these two solutions so that at x = the wavefunctions 
vanish. A linear combination for which both spinor components vanish at x = 0, i.e. au\ + Pu 2 = oti)\ + (3v 2 = 
yields vi/ui = V2/U2, hence an equation for the eigenvalues 



ii+ v /|A(+/c,fc y )| 2 -e 2 _ A(+k,k y ) 



-ii + ^/\A(-k,ky)\ 2 -i 2 A(~k, ky) 

Its solutions are readily seen to be e = —sign(k y )A(k 2 — k 2 )/k F . In terms of the incidence angle £, k y = k F sm(, 
k = k F cos £, the eigenvalues are 

= — sign(C)A cos(2C) vpA v sm£ (9) 



FIG. 1: Surface (vertical line) and angle £ where the spectrum is e = (full lines); in the absence of A y these lines would be 
at £ = ±7r/4 (dashed lines). The range for which e > is shown as the hatched area; the spectra spans the range e — up to 
e = A as shown 



Note that the spectrum is not symmetric in k y or in £ (it is in fact antisymmetric) resulting in a finite surface current. 
Fig. 1 shows the angle £ where ef = (full lines) and the range for which > 0. The velocities de^/dky are 
positive for both ±k y branches, i.e. the surface states are chiral. This property leads to quantization of Hall effects, 
as discussed in section IV. 

We note that self consistency would imply that A' = at x = the eigenfunctions would then be 
~ exp[— L A'(x')dx'\ sm£|/«ir], resulting in a very similar dependence on Note also that quasiparticles in the 
bulk have a spectral gap |A(fc, k y )\ which for any given £ is higher than the surface states (neglecting the Doppler 
shift). Impurities, however, may destroy k y conservation and scatter high energy (> A') surface states into degenerate 
bulk states. When impurity scattering is essential (e.g. section IV) our results apply only when these excitations can 
be neglected, e.g. at T < A'. This restriction is not needed at the (110) surface where the whole surface spectrum is 
< A', i.e. below the lowest bulk state. 

The decay length of the surface states becomes, using Eqs. £019, (3/T 1 = [{m/k)y/\A\ 2 - e 2 ]- 1 = sinCI with 
£' = vf/A 1 . Since \ui\ — \vi\,i — 1,2 (Eq. [JJ the normalized eigenfunctions are 



where L y is the length of the surface. It is remarkable that |iif(r)| = |«c( r )l for all £, i.e. for all energies of the surface 
states, implying maximal electron-hole mixing. As noted above, a (110) surface has the same solution (|l()(l with £' 
replaced by £ = vf/A. 

We note that in general the spinor Eq. @ can be decomposed in terms of eigenoperators rj^, rj^ i where 




sign(k y )u(r) 



(10) 




(11) 
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leading to the diagonal Hamiltonian 

H = Y,I dxe^l^+rtl^a-lWHx^ + Hx)] 2 ] (12) 
c J 

with being x dependent via the Doppler shift. The spectrum has exact particle- hole symmetry, i.e. for each 
eigenvector u, v with eigenvalue e there is an eigenvector — v*, u* with eigenvector — e. The form Eq. (|12[l incorporates, 
however, both ±e states and its sum is therefore restricted to > 0. 

We consider next a d + is state at a (110) surface with an order parameter 

A(p x ,p y ) = Ap x p y /kp + iA s . (13) 

Eq. © has then the solution e = sign(k y )A s , i.e. 

e = signCA s vpA y sin £ . (14) 

Positive eigenvalues are now at k y > (for weak Doppler effect -Ui?|^4j,| < A s ) with a weak dispersion due to the 
Doppler term. Note in particular that the spectrum has a gap, i.e. no e = states; hence to probe these states one 
needs either high voltage or high temperature T > A s . This d + is state corresponds to a (110) surface at which it 
breaks both parity and time reversal. At a (100) surface the state d + is state is symmetric under reflection and in 
fact has no surface bound states. Hence tunnelling data at the (100) may distinguish between d + id' and d-\-is states, 
i.e. the d + id' state shows a weak structure at a bias w A while a d + is state has no effect at all. The magnetization 
data 3 shows an effect for both (110) and (100) surfaces, supporting a d+ id' state for YBCO. 



III. SPONTANEOUS MAGNETIZATION 

The d + id' or d + is order parameters break both time reversal invariance and reflection along the surface, hence 
they allows surface currents (d + is refers to (110) only). The current density parallel to a surface (the y direction) 
and the charge density are, 

n edge = ~£<*i(r)* s (r)) = f £k(r)| 2 (15) 

C 

where d is the interlayer spacing, (7^ s 7^, s ) = [1 + cxp(e^/T)] _1 and |uf (r)| = |«c( r )l were used. The expression for 
jedge can also be obtained from Eq. l|12f) by j e d ge = c8H / 8A y {x). In addition to the explicit T dependence in l|15|) the 
order parameters are T dependent as Aw A \/t , A' w A' q ^/t where t — (T c — T)/T c ; hence £ w foA/^i£' ~ Co/v^"- 

In principle the current has also a diamagnetic term (e/c)n e d ge (x)A y (x); the ratio of this term to the London term 
{cj Air\ 2 )A y (x) is l/[(kp£' Q ) 2 T] where r = (T c — T)/T c . Hence the effect of this diamagnetic current is small except 
very near T c , i.e. for r > (^f^o) 2 ( or T > (^f£o)~ 2 for the (110) surface). In the range where the order parameter 
fluctuations exceed its mean value mean field breaks down; this range, which is between r < l/kp^Q m 2-dimensions 
and r < (l/fcp^o) 4 m 3-dimensions is excluded in our analysis. 

We consider first d + id'; the factor k y tanh(e^/2T) is symmetric in k y , therefore within the integration in Eq. 
(fT5|l on the > range (Fig. 1) the ( < segment can be shifted into a ( > one so that a complete (0, n/2) 
range results. In terms of the density n = k F /2nd and the T = penetration length Ao (A w \ /y/r) where 
Aq~ 2 = iirne 2 / mc 2 — 2kpe 2 /mc 2 d we obtain 

/ ~.\ _ 2 ^o r' 2 ^ „„„ a „:.a ^-2* sm ae r A cos 2 C + ( e / c )^ sin CA, ( s) 1 

2T 



— 3ed B e{x) = 4^2 / dCcosCsm 2 Ce- 2a:sinC/ «'tanh[ 



A Jo 

ek 



n edge (x) = ^ / V2 dCcosCsinCe- 2a:sin ^' (16) 

where the rapid oscillatory sin 2 kx is replaced by its average i . Note that for A = or A' = all angles £ are allowed 
in the solution of Eq. (JSJ) and then the current vanishes. This demonstrates that BTRS leads to current carrying 
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surface states. We note also that if the Doppler shift ~ A y (x) is ignored the integrated current j e dge(x) vanishes, 
unlike the p wave casein. 

The response of the condensate to j e dge involves the London terms as well as coupling to the scalar potential at 
the surface; the latter terms are small as at low T or vanish at T — > T c (see Appendix C). London's equation 

with jedge (%) as a source term is then, 



- V 2 A y ( X ) = {-(l/\ 2 )A y (x) + (4TT/c)je d ge(x)}6(x) 



(17) 



where 6{x) is a step function. This assumes a thick film, i.e. no dependence on the z direction; the thin film limit is 
considered below. For a thick film the condition of no external field at x — > — oo implies H y (Q) — 0. Eq. I|17|l is then 
solved by the Greens' function 



G(x,x') = -(X/2)[exp(-\x-x'\/X) + exp(-\x + x'\/\)] 



(18) 



which satisfies the boundary condition d x G(x,x')\ x —o — equivalent to H y {x = 0) = d x A y \ x= o = 0. A y {x) then 
satisfies an integral equation 



Ay(x) 



00 A 

<'A§ Jo 



r/2 



d( cos £ sin ( 



dx' tanh[ 



A cos 2£ + (e/c) vp sin ( A y (x' 



2T 



)][ e -|x-a'|/A + 



e ~\x+x'\/\-i e -2x' sinC/f' 



(19) 



The doppler shift, as shown below, is significant only very near T c or at very low temperatures. Neglecting first the 
Doppler shift and at T — > T c Eq. l(T^|l becomes 



,7r/2 

A y (0) = (2(/. AA/7rA^r c ) / rfCcosCsin 2 Ccos2C(2sinC + £'/\Y 

Jo 



(20) 



The total spontaneous flux is $ = A y (0)L y where L y is the length of the boundary. We consider 20o_L l) AA/7rA o ! T , c as 
a flux unit, e.g. for L y = 2cm& and typical YBCO parameters it is « lO 5 0o. This flux unit is weakly temperature 
dependent since AA w AoAo is finite at T — > T c . The ratio $ = — ^ / (2<j)[)L y \/S./iT\ 2 ] T c ) is plotted in Fig. 2; it varies 
between A/12£' at A < £' (weak BTRS) to C'/12A at A > £' (strong BTRS) with a maximum of 0.014 at A w For 
a (110) surface replacing ^' by ^ (considering only £ -c A) we obtain $ = ^A'/12AA, much smaller than for a (100) 
surface. The reason for the dominance of the (100) surface is the steeper spectra e ~ A for this case. The result that 
$ is weakly temperature dependent at T — > T c is consistent with the spontaneous magnetization datai; more details 
on the data follow in section VI. 



0.015 




FIG. 2: Spontaneous flux for a (100) boundary in thick films (£' < d) in units of 2(f>oL y \A/iT\oT c 
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TABLE I: Spontaneous flux <E> [flux in units of — 0o-L y AA/(7rAoT c )] for various surfaces, weak (£' 3> A) or strong BTRS (£' <C A) 
and various temperature limits. Comments: (i) If thin film is not specified the entry corresponds to thick films with thickness 
d 3> A,f;' [(100) surface] or d ^> A,£ [(110) surface], (ii) All entries correspond to d + id' except the d + is one which refers only 
to (110) surface; for thin films it is the same as (110) d + id' (except a factor -5 in the T — > T c line), (iii) T — + T c entries for <E> 
exclude paramagnetic anomaly regions which are given in the last column, (iv) BTRS which sets in at a temperature T' c <C T c 
has <3> values corresponding to modified temperature intervals. The only paramagnetic anomalies in this case are (100) T < T s 
and (110) thin film, (v) The fluctuation region (e.g. r < l/kp£o m 2-dimensions) is excluded, hence the region r = (T c — T)/T c 
in the last column is relevant only if it is a larger one. 



geometry 


temperature 


case T" <C T c 


l> 


p£Lr£LHi3,gii6tic cLnom&ly 


firing c' ^ \ 
(,iuuj t, s> a 


1 — > ±c 




A 
15?' 






T < T c 


T < Tc 


( 2 V2) 6 ?'A 




(100) f < A 


T^T C 




Jl 

12A 


\ X ) 




T< T c 


T < T' c 






(110) 


T^T C 


T -> 


e 2 

12A£' 






T s < T < T c , A' 


T s < T < T' c 








T < T s 


T <T S 


An 


T as £-T 


(100) thin film 


T^T C 




A 


'-(*)" 




T< T c 


T < Tc 






(110) thin film 


T^T C 


T -4^ 


15C'A 


Doppler dominated 




T< T c 


T< 2^ 


i 2AT C 
3?A 


Doppler dominated 


d + is 


T^T C 


T -> 








T< T C ,A S 




T c 
A 


ls ~ A Jc 



At low temperatures T <§; T c the result for the (100) surface is of the same order as that near T c while for the (110) 
$ is enhanced upon cooling, becoming at T <C T c , A' of the order of $ w £/A. The various limiting forms of $ are 
collected in table I. 

We consider next the results with the Doppler shift. For the (100) surface and £' ^> A the kernel G(x, x') is localized 
at x ~ x' so that A y {x') can be replaced by A y {x) in Eq. . Near T c we expand the tanh and obtain a term which 
modifies 1/A 2 , i.e. 

where j y °\x) is the current in the absence of the Doppler term. Very neat T c , the effective London length X e ff where 
1/Agif = 1/A 2 — A'/2TAq becomes imaginary so that there is no Meissner effect, i.e. a magnetic field can penetrate 
into the bulk. Hence a sharp sign change of A y (Q) from paramagnetic to diamagnetic is expected at r « (A'/A) 2 . 

For T <C T c we obtain from ffl that ^UpA^O) « A' -C A, i.e. the Doppler shift is negligible. 

For the (100) surface and £' -C A the a;' integration in Eq. (|19fl is limited to hence we can replace by 
Aj,(0) to yield 



Ay(x) = % <cosCsinCtanh[ ?!— 1-£ LJtii e -*/\ (22) 

7TA ./n 2J 



At T — > T r this becomes 



d-S§|K( ) = 4 0) (23) 

hence the response changes sign at r < (£/A) 2 . For T <C T c the Doppler term can be neglected l^-A^O) w 
AoCC'/^ 2 ^ Ao- For the (110) surface the form (|23|l applies with £' — > £ (considering always £ <C A) which does not 
affect the left hand side of Eq. hence a paramagnetic anomaly at r w (£/A) 2 . 
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FIG. 3: Spontaneous flux at a (110) surface for A' = (dotted lines) showing a critical temperature T a . For A'/A = 0.01 it 
shows enhancement below T s . 



A remarkable feature of Eq. (1221) is that it allows spontaneous magnetization for the (110) surface even if A' = 0, 
as studied earlie»ii*i£. The critical temperature can be deduced from Eq. (|22[) (with A — > A') by assuming a small 
probing A' and look for the A y response, which from Eq. (|23|l diverges at T s — vf/3Aq w (£/A)T c <C T c . Furthermore, 
at T = Eq. Q2) yields 



Ay (0) = -rL si gn{A y (0)} 

7TAo 



(24) 



hence a spontaneous magnetization flux of ±0o£ a /7rAo. 

In Fig. 3 we show the low temperature form of ^4 y (0) for (110). For A' = it shows a spontaneous magnetization 
(dotted lines) below a critical temperature, while for A' ^ it shows enhancement near T s where it joins one of the 
low T branches. In comparison the (100) flux depends weakly on temperature and is much stronger than that of (110) 
at least at high temperatures. 

We consider next the thin film case, for which London's equation is 



X7 2 A(x, z) = d[-^A y {x) + —j(x)]9(x)S(z) 
A z c 



(25) 



where d is the film thickness. Assuming that one can Fourier transform A(x,z) into A(q,k), integration of the k 
dependence ~ (q 2 + k 2 )^ 1 yields for A y (x) — A y (x, z — 0) 



A y (x) 



dq 
2^\q] 



dx'e^-^dl-^Ayix) + —j(x)]6(x) . 
A c 



The q integration then yields 



Ay(x)-Ay(0) 



dx' In 



- 1 47T 

d -Ta A/ * + — J'(a0 
A z c 



(26) 



(27) 



implying a slow decay of A y (x). While a solution for A y (x) appears difficult to obtain, the value of A y (0) is readily 
noticed from the boundary condition. The absence of external filed requires a finite H y (x = 0) for the thin film 
geometry. Hence to avoid divergence of dA/dx\ x= o where 



dA(x) 
dx 



dx' 7r 1 . , . 4-7T 

,- sl =/ —d[—Ay(x) j{x)\ 

x A 



(28) 
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FIG. 4: Doppler induced paramagnetic anomaly near T c for thin films with a (100) surface and A'/A = 0.1 



one must have A v (0) = \ ^rj(0), i.e. 

Ay (0) = / dC' cos C sin 2 C tanh 5=-^ (29) 

which interestingly has the same form as Eq. 1)221) except that here it is valid for all £'. In particular, when the 
Doppler shift can be neglected we obtain $ = A/15£' which in Fig. 2 is the tangent line to the thick film curve at 
the origin. Hence we can define two regimes: Weak BTRS with A/£' < 1 where the spontaneous flux is T and d 
independent, and strong BTRS with A/£' > 1 where film thickness matters, with the thin film showing a stronger 
effect. For strong BTRS T dependence is induced as £ ' < d changes to the thin film case £' > d as T — > T c . 

Consider now the Doppler shift for thin films; expansion near T c yields 1/A e // = 1/A 2 — A'/2TAq which as above 
changes sign at r w (£/£') 2 j i-e. a paramagnetic anomaly. This temperature is the same as for the thick film case 
except that here it is valid also for £' < A. Hence for A'/A = 0.1 we can have an anomaly at an accessible temperature 
of (T c - T)/T c w 1CT 2 , as shown in Fig. 4. For T < T c the Doppler effect is small. 

For a (110) surface the scales of A y (x) are A, £, hence the the thin film situation applies when d <C £, which is 
more difficult to achieve. Near T c we obtain A y (0) = j^A'A > which is paramagnetic, while at T C T c we have 

Ay(0) — igj^rj'. A spontaneous flux even with A' = is possible also here as in the (110) thick film case. 

All the various forms for the magnetization and Doppler effects are summarized in table I. The table also considers 
the possibility that BTRS sets in at a temperature T' c <C T c . In this case £' diverges at T' c so that near we have 
£' ^> A while at T <C T' c also £' < A is possible, resulting in a temperature dependent <I> for this strong BTRS case. 

Finally we consider the d + is case. Here only (110) is relevant and £ <C A for high T c materials. The edge current 

is 

-ie dfl e(x) = % r / ^CcosCsm 2 Ce- 2 - i "^ tanh[ A, + ( e /c)^sm C A,(,) l 

Near T c Eq. 1)20(1 applies with Acos2£ replaced by A s and £' — * £, hence $ = A S /2A = £/2£ s . At low temperatures 
the tanh is replaced by 1, leading to $ = T c /A. The various results are given in table I. 



IV. QUANTUM HALL EFFECTS - A SURFACE APPROACH 



In this section we study a surface formulation of quantum Hall effects (QHE). For the usual charge conduction, 
in the absence of an external field, and given that the surface fields decay in the bulk (Meissner effect as shown in 
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section III), Ampere's law implies a zero net current, i.e. a net Hall conductance <j xy = 0. We focus therefore on spin 
and thermal Hall effects. These were shown to be quantized first by a network model simulations^ and then by the 
relation to edge states^. The d + is state has no surface states near e = 0, hence no Hall effect within linear response; 
there may be a response when the voltage exceeds « 2A s /e. We consider therefore in this section only the d + id' 
case. 

The main ingredient is the chiral nature of the surface states. These states have two branches whose spectra vanishes 
at k y = Q = ±fci?/\/2- Linearizing near this point the spectrum for k y = q ± Q is 

e = vq (31) 

where v = V8A/k F for a (100) surface while v — \^8A'/kp for the (110) (up to a small Doppler shift; the actual 
value of v is not important for the eventual result for the Hall conductance). We wish to rewrite the surface modes 
in terms of two branches 771^,772,9 where 

m,q = Vq+Q] 
m,q = Vq-Q] 

m-q = -vI-qi (32) 

with the ± signs needed for continuity (see below). Thus, instead of two q > branches with two f , | spins we have 
now two branches, each of a single degree of freedom, with both q > and q < 0. The transformation of Eq. 
with the eigenfunctions i|I0|) is 

'*t( r )\ / 1 i\ / m. q e l(Q+q)v \ (\ -i\ f m, q e- l{Q - q)y \ 

+ \=J2^ x n )+f-q(x)[ } 

*t(r)/ t^o V-l lJ \72,-,e- i W + «)«/ \l 1 / V7i,V (Q " 9) V 



(33) 



where the two terms correspond to the two branches near ±Q and 



/,(*) = ^(^-(Q + ^^e-^+'l/^ . (34) 

Note e.g. that the Q ± q Fourier components of ^(r) are — ?7i l9 for both ±q so that the ± signs in Eq. (|32|l are 
needed for continuity of the Fourier transform. We can therefore define fields with continuous Fourier transforms 

*i(r) = xy< o+9) v*/ 9 o«o 

q 

* 2 (r) - Y. el{ ~ Q+q)V ^qf-q( x )- (35) 



Eq. I|33|l becomes 



or in terms of spinor 



* T (r) = ^i(r) + * 2 (r) 

*t(r) = -*i(r)+* 2 (r) (36) 



*W = (37) 



we have for the spinor Eq. (|f f \ 



.*a(r). 



*(r)= I 1 = ( ! |)* (r) - (38) 
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The kinetic energy has the form 

Hk = X><M^M + 4, q m, q ] = [ dxdy¥(r)[-ivd v Mr) . (39) 
q J 

A crucial observation for QHE is the role of impurities which in general has the form 

n imp = J dxdy^ (r) n ^(r)V(r) (40) 
Within the subspace of surface states this interaction becomes, using Eq. (|38() 

U lmp = j dxdy¥{T)T^{r)V{r) (41) 
This impurity potential can be " gauged" awayi^ by a transformation 

*(r) - exp[(i/r;) C V(x, y')dy'n]$(T) (42) 



which eliminates the impurity potential. Hence transport of chiral states is equivalent to that of a pure system. 
Chirality implies no channel for backscattering, hence impurities are indeed expected to be ineffective. As noted in 
section III, (100) surface states with energy > A' may mix with bulk states by impurities. The QHE is then limited 
to temperatures T < A'. 

To evaluate the spin Hall conductance, we define a spin voltage V s (x) (x is a coordinate perpendicular to the edge) 
such that (h/2)p s (x)V s (x) is the coupling energy density to a density p s (x) of h/2 spins. This can be represented by 
a Zeeman term with V s (x) = eB z (x)/mc where m is the electron mass; the corresponding force in the x direction is 
(e / mc)dB z (x) / dx . The unit of spin conduction, in analogy with e 2 /h of the charge conduction is (h/2) 2 /h = H/8tt. 
Hence the spin Hall conductance (J x p v m = I s /V s is 

a s / y m = 2 -^sign(AA') (43) 
The thermal Hall conduction is derived in a similar way from the heat conduction of an ideal gas, yielding 

^¥ = ^ si § n ( AA ') ( 44 ) 

Hence K xy /T is also quantized in this weakly disordered system. 

We reconsider now the effect of disorder on the spin Hall conductance. Imagine many random Hall systems, each 
with their own localized chiral states which are weakly coupled. If the couplings are too weak we expect no currents 
between the systems so overall <J s x V y n — 0. As the coupling strength increases we expect a finite current to circulate 
around the ensemble of grains, leading to Eqs. (|43I44|I . The transition is in fact induced by disorder: For weak 
disorder the argument of Eq. (|42[1 holds and the Hall coefficients have their quantized values, Eqs. (|43I44|) . As 
disorder increases, opposite chiral channels get coupled leading to formation of localized chiral loops which eventually 
have an insulating behavior, i.e. the Hall conduction vanishes. This quantum Hall plateau transition^ has been 
simulated by a network model, showing a novel type of QHE criticality. 



V. QUANTUM HALL EFFECTS - A BULK APPROACH 



We consider next the effective action of a bulk d + id' superconductor and derive its (charge) Hall conductance 
a xy (q,uj). We assume a thin film situation with the scalar and vector potentials (j), A being z independent, as well as 
A z = 0. In terms of the Nambu spinors Eq. © the off-diagonal Hamiltonian J d 2 7"i/^(r)ft,A'0( r ) (Eq- IA6(1 is 

h A = -[A(-fi£ + d 2 y )n + A'd x d y T 2 ]/k 2 F (45) 

and we neglect terms with V# << fcj?. The action in presence of the electromagnetic potentials A, (p is then 

S = J d 2 rdt^(id t - r 3 e(p) -h A - T,)tp 

£ = T 3 (ao + a 2 /2m) + a ■ p/m - iV ■ a. /2m (46) 
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where e(p) = (p 2 — kp)/2m and we introduce the gauge invariant potentials a = ^'VO — eA and ao = ^-j^O—etp. The 9 
derivatives arise from the spinor transformation ^(r) — > exp[ir 3 9 (r)/2]^(r). Integrating out the fermions (Appendix 
B) and expansion to 2nd order in a, ao leads to the effective action 

d 2 qdbj 



s eff= I 7T7^-- p ^(q I w ) a M(q J w ))^(-q. w )) ( 47 ) 



At T — and q, u> — > we obtain Poo = Ao (density of states which is Ao = m/2n in two dimensions), Py = — Aoc 2 
where c s = vf/V2, while Poj(q) — isign(AA')eoijqi/ (4-k) and eoij is the antisymmetric unit tensor. The latter term 
reflects BTRS and is derived for A' < A. 

Integrating out the phase 9 we obtain the effective action in terms of the electromagnetic potentials A, ip 

a r* i 2 f d 2 qduj c 2 s q 2 r . , . l2 % , . . 

S eff {A,(p} = e -TTT-^-i -9Z2 2"[ Po °l^(q' w )l " l- e Oyft^(q,w)Aj(-q,-w) 

J (z7T) C^q z — 47T 

+0(c 2 |A| 2 )] - P 00 (^) 2 |A(q,c)| 2 } (48) 

The total electromagnetic action includes also the Maxwell terms Sm = J d 2 rdt(E 2 — H 2 )/8ir. A may also be 
integrated out, using V • A = and A z — leading to the effective action 

X(g,w) c s9 ~ w o7re z (»7rc)^ (c^ — lo^Y 

d c s w2( ? 2 ^ 2 dq 2 

X(q,w) = Poo— + 75 x 2 2 p ^ 2 2 2T + q — S - ? ~ o — 2 ' 49 J 

<r (87r) z c' 2 Poo(cfq z — ur) onc^e* S^e- 2 

The coefficient of \4>\ 2 vanishes when g — ► at the plasma frequency uj p = {Aime 2 /m) 1 / 2 = c/Ao; there are no acoustic 
plasmons. 

The Hall current J y is identified by a functional derivative with respect to A y leading to the Hall coefficient 

a xy (q,w) = sign(AA')-f- 2 j q2 2 (50) 

Transport is defined by taking the q — > limit first, i.e. o" xy = 0. Hence the conventional Hall coefficient vanishes, as 
expected from Galilean invariance 16 . A limit in which uj — > is taken first yields a quantized "static" conductance 
e 2 /2h which was argued to correspond to a xy 7^ in presence of a boundary^. In absence of an external magnetic 
field, and given a spontaneous magnetization decaying in the bulk (as shown in section III), Ampere's law yields zero 
total current, hence a xy = 0; this is valid also with a boundary and an external electric field. It is intriguing, however, 
that a xy {(\, uj) has a nontrivial structure and space resolved measurement of a Hall current could then probe the full 
Eq. H5U[1 . We note that a result similar to Eq. (|50|) was obtained for superfluid 3 Hei£. 



VI. CONCLUSIONS 



We consider now in more detail the experimental data on the spontaneous magnetization^. The data shows that 
for a YBCO disc with a perimeter of L y w 2cm the spontaneous magnetization is temperature independent in the 
range 80-89K and is also thickness independent in the range 30-300nm with a value of «370o- Taking AA w A0A0, 
their T = value, and typical YBCO parameters we find $ w 10 -3 . The temperature and thickness independence 
indicate weak BTRS with £' > A. For either thick or thin films we estimate A/£' w 10 -2 or A'/A w 10~ 4 . We 
propose therefore that increasing the ratio A'/A, e.g. by using overdoped YBCO^, one can enhance the spontaneous 
magnetization up to a maximum of ~ 10 3 </>o when A'/A rs 0.01 within the thick film regime. 

For strong BTRS, A/£' > 1, the film thickness matters, i.e. we expect a temperature dependence due to the 
crossover from thick to thin film regimes at d w f as T — > T c . For thin films (d < f < A) we obtain $ = A/12£', i.e. 
for YBCO the total flux can reach 10 5 A'/A0o per cm of boundary, much higher than thick film values. The situation 
of a strong BTRS with thin films is interesting also as being the most likely one to show the paramagnetic anomaly 
(Fig. 4) at a temperature « T c [l - (£/£') 2 ]. 

In conclusion, we have shown that surface states of a d + id' superconductor lead to spontaneous magnetization 
which is T independent and thickness independent for weak BTRS, A/£' < 1, in accord with the data-. For strong 



12 



BTRS with A/£' > 1, as expected in overdoped YBCO^ , a crossover from thick to thin film behavior can lead to T 
and thickness dependence, as well as to an observable paramagnetic anomaly near T c . We have shown gapless chiral 
surface states for the d + id' state which lead to quantization of the spin and thermal Hall conductances. The d + is 
state has surface currents only at the (110) surface; its surface excitations have a gap and therefore insulating; i.e. no 
nontrivial quantization of Hall conductances. For the charge Hall conductance we find a vanishing transport value, 
however the structure of cr XJ/ (q, u>) has an unusual form which exhibits the Goldstone mode of the superconductor. 
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APPENDIX A: HAMILTONIAN FOR D WAVE SUPERCONDUCTOR 

We derive here the interaction term for a d+id' superconductor. A general interaction Hamiltonian in terms of a 
Nambu spinor Eq. J2J is 

Hint = - J tft( ri ) T 3^( ri )tft( r2 ) T3 ^( r2 ) F ( ri - r 2 ) (Al) 
The order parameter has the form 

e ^3 ri A(r l5 r 2 ) = (*(n)*t( r2 ))v(n - r 2 ) (A2) 

where the phase 8(r) depends only on the center of mass coordinate r = (r x + r 2 )/2. The factor A(r 1; r 2 ) may be 
complex, however, its real and imaginary components are determined by the interactions and their ratio is not allowed 
to vary in space. A d wave superconductor is defined by a momentum dependence k 2 — k 2 for the relative coordinate, 
i.e. 

A(n,r a ) = A(r) J ^~^{kl - k 2 ) = A(r)(-df + d 2 )S 2 (p)/k 2 F (A3) 
where the relative coordinate is p — Y\ — r 2 = (£, rf). The mean field Hamiltonian is then 

H&f = / d 2 rd 2 p^(r + l p ) e ^r^(r - ^p)A(r)[9f - d 2 }S 2 (r)/k 2 F (A4) 
After partial integrations, 

= \J d 2 r^{v)[{dl - d 2 )A(r)e ie ^]n^(r)/k 2 F 

- J d 2 r^(r)[d x A(r)e lf>ir)T *d x - d y A(r)e i6 ^ Ts dy]ri^(r)/k 2 F (A5) 

where in the first term (d 2 — d 2 ) operates only within the [] brackets. A d' component corresponds to A' (r)d^d v 6 2 (p) 
and similar analysis can be followed. We assume here that all gradients are small, i.e. |V6>|, |VA|/A <C kp, hence 
with the transformation ^(r) — > exp[iT^9(r) /2\^(v) yields the off-diagonal Hamiltonian J d 2 r^ (r)hA'^(v) where 

h A = -[A(-d 2 x + tf)n + A'd x d y T 2 ]/k 2 F (A6) 
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The issue of gauge invariance is of some interest. The full interaction form Eq. (|A1|) is manifestly invariant under 
^(r) — > exp[ir 3 J T A(r') ■ <fr']\l/(r). Wether the mean field form is also gauge invariant is a matter of some debat o 21 ' 22 . 
From the definition Eq. (|A3(I it seems that 

/ri pr 2 
A(r') • dr']exp[iT 3 / A(r') ■ dr']A(ri, r 2 ) (A7) 

and then Ttf^f is gauge invariant without having explicit A(r) dependent terms. This, however, implies that the d x , d v 
terms do not follow the usual substitution law as in the kinetic term. For the present work this issue is irrelevant 
since we neglect these terms altogether, i.e. we assume |V0|, |VA| <C /cf- 



APPENDIX B: DERIVATION OF P i} 

We derive here an effective action for a d + id' superconductor in terms of the gauge invariant potentials a^{q, ui), 
Eq. I|47|l. Integrating out the fermionic variables in the partition function we arrive to the following action: 

Z = J D$e lS (Bl) 

S(<f>) = -iTrlnG^ 1 , G-^Gq 1 -^ 
Gq 1 = id t - r 3 e(p) - h A 

We are interested in the long wavelength limit; also the order parameter is taken at the extremum of the effective 
action with only phase fluctuations. We retain the first and the second order in E to derive an expansion of the effective 
action in the fluctuating fields a^q^uj). The expansion corresponds to a one loop calculation with the coefficients 
Pfj, M in Eq. (|47|l given by (latin indices stand for space coordinates) 

P 00 (q W ) = l -T ]T Tr[G(p,u')r 3 G(p + q,u' + u)r 3 ] (B2) 

p.LU' 

PijW) = ^T^2Tr[G(p,LJ , )G(p + q,u , + w)p i (p + q/2) j ]-^5 u (B3) 
Pojiw) = ^TY,Tr[G(p,u')TzG(3p + q,J +u) Pj ] (B4) 

P,LJ' 

The diagonal time polarization operator Poo depends weakly on temperature and therefore in the limit of small 
momentum and frequency q — > 0, uj — *■ is given by its T — value, i.e. the mean-field compressibility, Poo(q) — 
No . The space components Pjj include the diamagnetic term and paramagnetic current correlator. In the limit 
q — > O.uj — > they give the mean-field superfluid stiffness; at T — > Pij(q) = —N c^ where c s = vp/y/2. Of special 
significance is the off-diagonal polarization bubble Poj which is responsible for the Hall effect. It is a topological effect 
depending (at least at small values of the d x , y order parameter) only on the sign of A'. In the same long wave-length 
limit we have Paj(q) = zsz<?n(AA')eoijQi/(47r) where eoy is the antisymmetric unit tensor. 



APPENDIX C: EFFECTIVE ACTION WITH BOUNDARY 



We study here the Hall term with boundary and show that its effect on London's equation is small at either T = 
or T — > T c . The electromagnetic response to the surface charge and currents couples in general the vector and scalar 
potentials A, cp with the Hall coefficient. We estimate this effect first at T = 0. The Hall term relates the current 
along the surface (jj direction) and the electric field dtp/dx in the x direction, i.e. 

1 d 2 47T dip 4-7T 

(j£ - ft?)4r(«) - Tc °*v-^ = --**.(*) (CI) 



The equation for (p(x) involves the Debye screening length = 1 /V 4me 2 <C £,',£,, 

& 



(^2 _ Q^^ix) = ^n edge (x) (C2) 
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The Hall term is neglected here as we wish to estimate the lowest order effect. The solution with dip/dx — involves 
the Greens' function Eq. (|18fl : at x <C £' it has dip/dx ~ Ad/£' while at x ^> 

v{x ) = _ r 2 dC cos C sin 2 Ce~ 2a; sin c/ « ' (C3) 



2nd^ J 

The ratio of dip/dx and jedge terms in Eq. I|C1|) is then l/(327rfc^^') <C 1; for a (110) surface replace £' — > £. Since 
kp^' > 3> 1 is the criterion for excluding the order parameter fluctuations very near T c , we can neglect the Hall 
term in London's equation IjClfl . 

We consider next T — > T c . The polarization function P 00 (B2) is obtained by replacing N — ► N(T) where 



with = a/ e 2 (p) + g 2 + g\ where 



gi = A(r)cos2C 

9l = A'(r)sin2C. (C5) 



The polarization function (B3) defines the temperature dependent London penetration depth 1/\ 2 L {T). 

We consider in more detail the Chern Simon (or Hall) part of the action which is the product of scalar and vector 
potentials (Eq. B4). We consider a superconductor that occupies the half space (x > 0) where the order parameters 
A, A' may become functions of x. The Chern Simon part of the action can be written in configuration space in the 
form 

S c - S = e 2 Tj2 J dr[h(rr')a (ruj)a y (rcj) + b 2 (rr')a (rcj) W (C6) 

P 

&2 ^') - ^Eft^.(r',r)ft(r')^ (C8) 

_ 2[^(#) - Eth&)\ 
tp(jr) ~ EE'(E' 2 -E 2 ) 

where E = E(r),E' = E(r'). Similarly as for the infinite system we can integrate out the Goldstone mode 

S b {A,ri = e 2 J dr[W(r) - ^{j^f A 2 {v) +b l {v) V {v)A y {v) +b 2 {v) V {v)^^-] (CIO) 

In this equation we took the polarization function at zero frequency which is legitimate for finite system (the effective 
momentum deviates from zero). The coefficients b\ — bi(rr')\ r —, r i , bi = b2{rr')\ r ^ r i are given as 

4 7 re 2 6 1 = JL^ffaT) (Cll) 

47re 2 6 2 = ±-f 2 ( x ,T). (C12) 
ncd 

The function fi appears only for a system with boundary and depends on the space derivative of the order parameters, 
while /2 at T = is the same as for an infinite superconductor . We calculate these function for T = and for T — > T c , 

h(*,T)= ^Zn[A(A') (1+5)/2 ] z/T = (C13) 

= 0.11^£^[AA'] ifT^T c (C14) 

Mx,T)= 1 z/T = 

0.2l4# if T ^ Tc (C15) 
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with (5 = 2/(1 + A'/A) and A(x),A'(x) > is assumed. Thus, in the limit T — > T c we we can write A « ^o/\/t, 
6i(r) = 0.11(AA7T c 2 )£Zn[AA']/2W and 6 2 (r) = 0.21(AA'/T c 2 )/2/icd. 

The next step involved the equations for electromagnetic potentials - generalized London equations by the variations 
of the total action (including the Maxwell part) over these potentials. We consider here the half-plane geometry, i.e. 
the superconductor occupies the x > half- plane. The nonzero electromagnetic potentials A y , ip obey the equations 

, 1 d 2 , Aire 2 dA v Aire 2 A . . . 

( A| dx~> - — h2 ~dx- ~ -dT blAy = *™«~<?) (C16) 

A—-„2 f), - A-rs>X 



1 d 2 w Aire 2 , dip i Aire 2 , t Q L N 4?r 

c 



where n e d ge , jedge are the edge charge and edge current densities (Eq. 1161). 

Using the expressions above for 61,62 we find that the Chcrn-Simon term affects the spontaneous magnetization, 
leading to an additional flux ~ (A'/A)t which vanishes at T — > T c , i.e. it is negligible compared with the other terms 
in Eq. I|C17|I which lead to constant magnetization as T — * T c , as shown in section III. 
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